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ABSTRACT 

High resolution far-ultraviolet spectra of the O-type subdwarf BD +28° 4211 were obtained 
with the Far Ultraviolet Spectroscopic Explorer to measure the interstellar deuterium, nitrogen, 
and oxygen abundances in this direction. The interstellar D I transitions are analyzed down to 
Ly i at 920.7 A. The star was observed several times at different target offsets in the direction 
of spectral dispersion. The aligned and coadded spectra have high signal-to-noise ratios (S/N 
= 50 - 100). D I , N I, and 0 I transitions were analyzed with cur ve-of- growth and profile fitting 
techniques. A model of interstellar molecular hydrogen on the line of sight was derived from H2 
lines in the FUSE spectra and used to help analyze some features where blending with H2 was 
significant. The H I column density was determined from high resolution HST / STIS spectra 
of Lya to be log N ( H I) = 19.846 ± 0.035 (2a), which is higher than is typical for sight lines 
in the local ISM studied for D/H. We found that D/H - (1.39 ± 0.21) x 10' 5 (2a) and O/H 

— (2.37 ± 0.55) x 10 -4 (2a). O/H toward BD +28° 4211 appears to be significantly below the 
mean O/H ratio for the ISM and the Local Bubble. 
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1. INTRODUCTION 

The abundance ratio of atomic deuterium to 
hydrogen (D/H) is a key diagnostic of light el- 
ement production and the baryon-to-photon ra- 
tio (77) in big bang nucleosynthesis (BBN, Walker 
et al. 1991). Significant efforts in recent years 
to measure D/H in high redshift intergalactic gas 
have generally yielded values in the range D/H 
= (2.5 — 4.0) x 10 -5 (O’Meara et al. 2001 and 
references therein). 

In galaxies, D/H is an important indicator 
of chemical evolution (Audouze & Tinsley 1974; 
Boesgard & Steigman 1985; Tosi et al. 1998), as 
D is easily destroyed in stars ( 2 H + p -+> 3 He + 7) 
at temperatures of only a few million degrees. In 
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fact, there are no known appreciable sources of D 
other than BBN. The interstellar medium (ISM), 
as the repository of stellar mass loss, holds the 
keys to unravelling the chemical history of D over 
the age of the galaxy. Measurements of atomic 
D/H in the ISM were first made by the Coperni- 
cus satellite toward bright OB stars within ~ 500 
pc of the sun (Rogerson & York 1973; see Vidal- 
Madjar & Gry 1984 for a review). Significant 
progress had been made using the Hubble Space 
Telescope ( HST ) to obtain some precise D/H mea- 
surements in the local ISM (Linsky et al. 1998; 
Vidal-Madjar et al. 1998), but these have nec- 
essarily been limited to very low column density 
lines of sight (N( H I) ~ 10 18 cm -2 ) where Lya can 
be used to measure the D I column density. 

The D/H ratio was measured at high spec- 
tral resolution (A/AA ~ 80,000) with the Inter- 
stellar Medium Absorption Profile Spectrograph 
(IMAPS, Jenkins et al. 1988, 1999) toward three 
OB stars ((SOri A, Jenkins et al. 1999; j 2 Vel 
and Q Pup, Sonneborn et al. 2000) first studied 
by Copernicus . These studies showed unequiv- 
ocally that there is a large spread in D/H (fac- 
tor of 3) on sight lines toward early type stars at 
distances of 300-500 pc. The IMAPS instrument 
was limited to observing very bright OB stars 
by its instrument design. At that time IMAPS 
could only observe the lower Lyman series through 
Lye A937 because its gratings were coated with 
lithium-fluoride (LiF). 

FUSE provides a new capability to study 
the entire Lyman series (except Lya) in much 
fainter stars than were possible with Copernicus 
or IMAPS, and hence D/H in the ISM within 200 
pc can now be studied systematically for the first 
time. This is important because the lines of sight 
to nearby stars are usually simpler with less blend- 
ing of different velocity components. FUSE has 
sufficient sensitivity to perform absorption line 
studies of white dwarfs as well as extragalactic 
objects. Both categories of targets were inaccessi- 
ble to Copernicus or IMAPS. The study of D/H 
throughout the Galaxy is one of the primary ob- 
jectives of the FUSE mission (Moos et al. 2000). 

In this paper we report the analysis of D, 
N, and O abundances toward the O subdwarf 
BD +28° 4211 using high signal-to-noise (S/N) ra- 
tio spectra obtained with FUSE. BD +28° 4211 
has a very high temperature (~82,000 K), lies at 


a distance of ~ 100 pc, and has a H I column 
density, IV (H I) , of ~ 10 20 cm" 2 . This column 
density is much higher than those of the other 
LISM lines of sight studied previously for D/H, 
and very similar to that of the D/H sight lines 
studied by IMAPS. In this 7V(H I) regime the Ly- 
man series shortward of Ly 5 A949 is required to 
study N(D I) . 

Because a measurement of D/H requires both 
accurate iV(D I) and N( H I) , we use FUSE spec- 
tra to measure total column densities for iV(D I) , 
7V(N I), and N{0 l)and HST j STIS spectra of 
Lya to measure 7V(H I) . Even in the regime 
N( H I) ~ 10 20 cm" 2 , the damping wings of Ly/3 
and higher Lyman series transitions are too weak 
to obtain an accurate 7V(H I) measurement. 

The observations and data processing are pre- 
sented in §2 and the stellar atmosphere model and 
synthetic spectrum in §3. The measurement of 
N( H I) is discussed in §4. In §5 we describe the 
measurement of N(D I) , 7V(N I) , and N( O I) by 
curve of growth and profile fitting techniques, as 
well as the analysis of interstellar molecular hy- 
drogen. This H 2 model is used in the D I , N I, 
and O I analysis. The paper concludes in §6 with 
a discussion of the results and a comparison with 
other recent D/H studies in the Galaxy. 

2. OBSERVATIONS 

BD +28° 4211 was observed by FUSE through 
the LWRS aperture (30"x30") on 2000 June 13 
for 2192 s (ObsID=M1080901), 2000 July 16 for 
16677 s (ObsID=M1040101), and 2000 September 
19 for 7864 s (ObsID=M1040105). One MDRS 
aperture (4"x20") observation was made on 2000 
July 17 for 24795 s (ObsID=M1040102). The first 
observation was a short exposure to establish the 
feasibility and safety of using BD +28° 4211 for the 
FUSE calibration program because its far-UV flux 
is very close to the bright limit of the FUSE de- 
tectors (F\ ~ 10" lo erg cm" 2 s -1 A" 1 ). The sub- 
sequent observations were made for the purpose of 
testing the acquisition of data with the Focal Plane 
Assembly (FPA) offset to different positions along 
the dispersion direction as a technique to reduce 
fixed-pattern noise when the observations are later 
coaligned and added. This procedure and the re- 
sults are discussed below. The FUSE instrument 
and its performance are described by Moos et al. 
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(2000) and Sahnow et al. (2000). 

The data were obtained in “histogram” mode 
because of the high far-UV flux of the star. In this 
data acquisition mode a two-dimensional spectral 
image is accumulated on board, eliminating any 
timing information within individual exposures. 
The observations were split into four intervals of 
approximately equal exposure time, each with a 
different X-offset of the FPA. The exposure time 
required for the last two LWRS FPA positions was 
obtained in part during observat ion M 1040101 and 
in part during observation M1040105. A shift in 
FPA X causes a corresponding shift in the location 
of the spectrum on the detector along the disper- 
sion direction. The spectral shifts and exposure 
times are given in Table 1. The FPA is fixed for 
the duration of each exposure so that there is no 
loss in spectral resolution relative to an observa- 
tion taken at a single FPA location. The offsets 
for the MDRS observations span a smaller range 
in FES X than those of the LWRS observations 
in order to provide sufficient margin in the FPA 
range of motion. This margin is needed so that 
MDRS peakups can properly adjust the FPA’s to 
correct for image motion caused b} small rotations 
of the primary mirrors on orbital timescales. For 
further information on this image motion see Sah- 
now et al. (2000). The shifts caused by this image 
motion, together with shifts caused by similar ro- 
tations of the gratings, result in a distribution of 
shifts clustered about the mean offsets given in Ta- 
ble 1. Even for MDRS, these spectral offsets span 
a range roughly twice the width of typical inter- 
stellar absorption lines or of the high-frequency 
components of the detector fixed pattern noise. 

Exposure durations varied from 430 s to 535 
s in order to minimize Doppler smearing of the 
spectrum by satellite orbital motion so that there 
is no degradation of spectral resolution. Most of 
the LWRS exposures were obtained during orbital 
night when 0 I and N I airglow emission is absent 
and that of H I is significantly reduced. When 
exposures having both orbital night and day seg- 
ments are considered, over 70% of the LWRS data 
was obtained during orbital nigh;. Over 82% of 
the MDRS data was obtained at night. 

The data were processed with the latest ver- 
sion of the FUSE calibration pipeline (CALFUSE v. 
1.8.7). For each channel, data for each exposure 
were aligned by cross correlation and combined on 


a common heliocentric wavelength scale. Part of 
the resulting spectrum for the LWRS observations 
is shown in Figure 1. The LWRS LiF co- added 
spectra have S/N~ 100 per 0.05 A resolution el- 
ement in the LiF channels and S/N~ 60 in the 
short-wavelength silicon-carbide (SiC) channels. 
The S/N in the MDRS spectra of BD 4-28° 4211 
are about 25% lower than for LWRS. These S/N 
levels are close to the photon-noise limits. The 
analyses presented in this paper uses primarily SiC 
channel spectra. Figure 2 shows small sections of 
the spectrum near several of the Lyman lines and 
demonstrates the high quality of our spectra. Our 
analysis of spectral features in the BD +28° 4211 
co-added spectra indicates that the spectral res- 
olution is consistent with A/AA ~ 20,000, the 
nominal spectral resolution obtained by the FUSE 
instrument (Sahnow et al. 2000). It is evident 
from Figures 1 and 2 that airglow has a negligi- 
ble effect on these spectra. Ly (3 is the strongest 
airglow transition in the FUSE band pass and 
the observed Ly/3 profile is only weakly affected. 
O I airglow, if present, would primarily effect the 
stronger transitions longward of ~ 950 A, but 
the O I column density is determined from much 
weaker lines that are unaffected. 

Figures 1-3 show that the saturated H I lines 
have non-zero residual flux (on the order of a few 
percent). The FUSE instrument has very low lev- 
els of scattered light (Moos et al. 2000, Sahnow 
et al. 2000), as shown by the zero flux level below 
the Lyman limit (see Fig. 1). The residual flux 
in the Lyman lines is probably the result of wings 
of the instrumental line spread function (LSF; see 
Kruk et al. 2001). The shape of the LSF and 
its potential effects on the line profile analysis is 
discussed further in §5.1. 

3. STELLAR MODEL FOR BD +28° 
4211 

MacRae, Fleischer, & Weston (1951) reported 
the first spectral analysis of BD +28° 4211 and 
classified it as an extremely blue star with pe- 
culiar spectrum (Op). The optical spectrum of 
BD 4-28° 4211 attracted their attention because 
of its very blue continuum and the presence of 
Balmer and He II lines. Greenstein (1952) esti- 
mated the distance to BD +28° 4211 and found 
out that the star was a subdwarf. BD +28° 4211 
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is now classified as sdO, based on the study of 
Moehler et al. (1990), in which the class is de- 
fined as having He II and strong Balmer absorp- 
tion lines. Although low resolution optical spec- 
tra show only these lines, high-resolution spectra 
reveal, in addition, many weak metal lines. Her- 
big (1999) obtained high-resolution optical spec- 
tra of BD +28° 4211 at the Keck I telescope and 
detected a significant number of narrow absorp- 
tion and emission lines (C IV, N IV, N V, 0 IV, 
O V, and Si IV), and also emission cores in Ha, 
H/?, and He II AA6560, 5411, 4685 lines. Us- 
ing the O IV absorption lines, Herbig (1999) set 
an upper limit on the star’s rotational velocity, 
usini < 4 kms -1 . Even though BD -1-28° 4211 is 
a well-known spectrophotometric standard star in 
the ultraviolet and optical (see, e.g., Bohlin 1986; 
Massey et al. 1988), its atmospheric parame- 
ters were investigated only recently by Napiwotzki 
(1993), Werner (1996), and Haas et al. (1996). 

Using H-f He NLTE stellar model atmospheres, 
Napiwotzki (1993) tried to fit the Balmer line pro- 
files in BD-b28°42irs optical spectrum with no 
success. He demonstrated that the effective tem- 
peratures obtained from fits of individual Balmer 
lines showed large discrepancies. This inconsis- 
tency, known as the Balmer line problem, was also 
observed in central stars of old planetary nebu- 
lae and DAO white dwarfs (see, e.g., Napiwotzki 
& Schonberner 1993; Napiwozki & Raunch 1994; 
Bergeron et al. 1994). To overcome this problem, 
Napiwotzki (1993) modeled the He I A5876 line, 
as it is a very sensitive diagnostic of effective tem- 
perature (T e ff), and obtained T e ff = 82,000 K and 
log(He/H) = —1.0. He also derived the gravity 
log g = 6.2 from the higher Balmer lines. 

By analyzing the optical spectra of BD 4-28° 4211 
and the central star of the planetary nebula 
S216, Werner (1996) concluded that the Balmer 
line problem was caused by the omission or in- 
adequate inclusion of metal opacities. For in- 
stance, he matched the Balmer lines (Ha-H<5) 
of BD +28° 4211 combining Napiwotzki’s (1993) 
parameters (T e ff= 82,000 K, log# = 6.2, and 
log(He/H) = —1.0) and the improved Stark- 
broadened CNO line profiles. He demonstrated 
that the effect of this added opacities increased 
the temperature in the deeper layers of the at- 
mosphere, and decreased the temperature in the 
superficial layers, therefore modifying the Balmer 


line formation. He noted that the low Balmer se- 
ries members cores are deeper as they are formed 
higher in the atmosphere, where the temperature 
is lower. On the other hand, the high Balmer 
series members were less affected because they 
originate from deeper layers of the atmosphere. 
This new temperature structure caused by the 
addition of more realistic metal line opacities re- 
solved the Balmer line discrepancy. Interestingly, 
Bergeron et al. (1993) had reached a similar con- 
clusion to that of Werner (1996) by analyzing the 
DAO white dwarf Feige 55 and using a LTE iron- 
blanketed model. 

For the purpose of this study, we adopt the stel- 
lar atmospheric parameters determined by Napi- 
wotzki (1993) as they are the best available. We 
compute a grid of H+He NLTE model atmo- 
spheres with the following atmospheric parame- 
ters: T eff = 78000 K, 82000 K, and 86000 K; 
log g — 6.0, 6.2, and 6.4; and log(He/H) = —0.6, 
— 1.0, and —1.4. We also compute a metal-line 
blanketed model including the following chemical 
composition: H, He, C, N, O, Si, S, Fe, and Ni. 
We then compute a grid of synthetic spectra from 
Lya to the Lyman limit using the TLUSTY pro- 
gram developed by Hubeny & Lanz (1995) and 
an upgraded version of the program SYNSPEC 
(I. Hubeny 2000, private communication), which 
incorporates the Lemke (1997) Stark broadening 
tables for hydrogen. We explored the model pa- 
rameter space {T e ff, logp, and log(He/H)) to eval- 
uate the uncertainties in the stellar H I Lyman + 
He II Balmer line profiles. §4 describes in more de- 
tails the use of stellar models to dermine the H I 
column density. 

The FUSE spectrum of BD +28° 4211 (Fig. 1) 
shows many stellar and interstellar lines. The stel- 
lar spectrum displays the Lyman series of hydro- 
gen from Ly/? up to Ly-9 and the He II Balmer 
series from 1084 A up to at least 942 A. The spec- 
trum also contains many stellar lines of highly 
ionized species, such as C IV, N V, O IV, O V, 
O VI, Si IV, P V, S IV, S V, S VI, Fe VI, Fe VII, 
Ni VI, and Ni VII. The strongest lines are the N IV 
A923.15 sextuplet, S VI AA933.38, 944.52 doublet, 
N IV A955.33, and O VI AA 1031.91, 1037.61 dou- 
blet. In order to identify potential blending of 
photospheric and ISM lines, we first performed an 
abundance analysis of BD -1-28° 4211’s atmosphere 
and then calculated a synthetic spectrum using the 
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derived abundances. Figure 3 illustrates the com- 
parison between the synthetic spectrum and the 
FUSE spectrum for the wavelength range 920-927 

A. 

BD +28° 4211 has several strong, narrow metal 
lines in the FUSE range and that made it pos- 
sible to register the wavelength scale of the stel- 
lar model with the observed spectrum. In par- 
ticular, S VI A944.523 is seen in both SiC chan- 
nels and has a deep, narrow core We used this 
line to align the model with the spectra. As ex- 
pected, the zero point wavelengths are different 
from channel to channel because it depends on the 
FPA position as well as the location of the star 
in the aperture. It varied between —60 km s -1 
and +60 km s“ x in our data, although larger off- 
sets were occasionally detected in the first year of 
FUSE operations (Sahnow et ah 2000). Once the 
wavelength offsets are corrected, the best estimate 
from the FUSE data of the velocity of the star rel- 
ative to the ISM is +34 ± 6 km s _1 , a value that is 
consistent with our more accurate determination 
of +31 ± 2 km s" 1 from HST spectra (see below). 

We increased the error to ±8 krn s _1 to take into 
account any remaining uncertainties in the FUSE 
wavelength calibration. 

4. DETERMINATION OF THE H I COLUMN 
DENSITY 

4.1. Hubble Space Telescope Observations 

The determination of the interstellar H I col- 
umn density along the sight line to BD +28° 4211 
uses an extensive series of archival observations 
taken with the Space Telescope Imaging Spectro- 
graph (STIS) on board the HST. The H I column 
density may in principle be determined using the 
FUSE observations of this star. However, most of 
the Lyman-series lines in the FUSE bandpass are 
on the flat part of the curve of growth. In addition, 
the damping wings of interstellar Ly/? and Lyy, 
the strongest H I transitions in the FUSE band- 
pass are significantly less prominent than those for 
Lya and comparable in strength with the stellar 
HI + He II profiles. Therefore, the only reli- 
able procedure for determining an accurate H I 
column density to BD +28° 4211 is by fitting the 
Lorentzian wings of the very strong Lya line. As 
illustrated below, the interstellar H I profile is 
much stronger that its stellar counterpart. 


BD +28° 4211 is a calibration source for STIS, 
and observations have been acquired over sev- 
eral years with the E140M echelle grating pro- 
jecting the two-dimensional spectrum onto the 
FUV MAMA (Kimble et al. 1998; Woodgate 
et al. 1998). All data were taken through the 
O'.' 2 x 072 aperture and cover the spectral range 
~ 1150 — 1750 A. This instrumental setup pro- 
vides a spectral resolution R = A/AA ~ 45,000 
(Av ~ 7 km s -1 FWHM), although the line spread 
function for data taken through this aperture has 
significant power in broad wings that extend well 
beyond this width (see Figure 13.87 in the Cycle 9 
STIS Instrument Handbook). We used eight such 
observations taken over the time period 1997 Sept. 
19 to 2000 Nov. 10. Exposure times for the indi- 
vidual observations were in the range 350-1800 s. 
Figure 4 shows a coaddition of the data sets used 
in our analysis in the region near Lya. This figure 
demonstrates the quality of the data and shows 
the importance of high spectral resolution for sep- 
arating the numerous stellar and interstellar lines 
in this region. 

4.2. Data Calibration and Reduction 

The E140M data have been reduced using the 
CALSTIS pipeline by two different approaches in 
an attempt to uncover sources of systematic error. 
One used the pipeline as developed by the STIS 
Instrument Definition Team, the other the stan- 
dard STSDAS pipeline distributed within IRAF 
(see Voit 1997). The CALSTIS pipeline removes 
an estimated dark level, applies a flat field, cor- 
rects small-scale distortions, identifies the spectral 
trace, removes scattered light, extracts the one- 
dimensional spectrum, and performs the wave- 
length and flux calibrations. We show below that 
both analysis approaches yield the same result. 

The scattered-light removal process is among 
the most important, particularly for our purposes 
where we are interested in the large-scale wings 
of the Lya profile. Because the strength of the 
echelle-scattered light can be significant in E140M 
spectra, the Bowers & Lindler algorithm (see 
Landsman & Bowers 1997; Bowers et al. 1998) 
was used to estimate and remove scattered light 
from the data in the pipeline procedures. We 
found that this algorithm does a good job of de- 
termining an appropriate zero-level, as evidenced 
by the small residual fluxes in the saturated cores 
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of strong lines (e.g., the Lya profile shown in Fig. 
4). 

After extracting the spectrum, several echelle 
spectral orders in the regions adjacent to Lya 
were combined using a weighted averaging scheme 
where the orders overlap. The flux calibration be- 
tween orders seems good for the E140M data sets 
used here. Even so, there can be some small fluc- 
tuations in the final co-added spectrum due to mis- 
matches at the edges of orders, although in general 
these can be masked out and do not significantly 
affect our analysis. 

4.3. Analysis Methodology 

The method for deriving the H I column den- 
sity along this sight line follows Jenkins (1971) and 
subsequent work (see discussion, particularly of er- 
ror analysis, in Vidal-Madjar et al. 1998; Howk, 
Savage, & Fabian 1999; Jenkins et al. 1999; Son- 
neborn et al. 2000). The distribution of optical 
depth, t(A), with wavelength for Lya is given by 
the product of the absorption cross-section, cr(A), 
and the column density. For H I at Lya this is: 

r(A) = a(X)N = 4.26 x 10“ 2O JV(A - A 0 ) -2 , (1) 

where A 0 is the Lya line center at the velocity 
centroid of interstellar hydrogen. The H I col- 
umn density is determined to be the one that best 
matches the above distribution of optical depth 
with the observed profile. This method only gives 
information on the total column of H I , since the 
details of the component structure are all within 
the strongly-saturated core of the Lya line which 
spans a velocity range of ±200 km s _1 . It is ex- 
tremely unlikely that the damping wings are af- 
fected by low column density H I absorption fea- 
tures. Any H I absorber that could modify the 
shape of the damping wings would have to have 
a heliocentric velocity \v\ > 400 km s _1 (see Fig. 
9 of Sonneborn et al. 2000). The broad damp- 
ing wings of Lya are not sensitive to the Doppler 
broadening parameter. We therefore have very lit- 
tle information about b from the Lya analysis. 

The derivation of 7V(H I) was performed with 
two independent, though similar, fitting proce- 
dures. Both vary the important free parameters, 
which for our purposes include A 0 and N ( H I) , to 
minimize the \ 2 statistic between the model, after 
convolution with the instrumental LSF, and the 


observed spectra (see, e.g., Jenkins et al. 1999). 
We note that the specifics of the LSF are rel- 
atively unimportant for N ( H I) as large as is 
found on this sight line (see, e.g., Appendix A 
of Howk et al. 1999). This insensitivity includes 
the possible presence of the broad wings of the 
STIS E140M LSF mentioned above, which are on 
a much smaller scale than that of the H I absorp- 
tion. 

We applied masks to exclude obvious stellar and 
interstellar absorption features from the fit. We 
also avoided fitting the very deepest portion of the 
Lya line core (effectively those regions with flux 
< 10% of the stellar continuum within ~ ±2 A 
of the line center) because of the possibility that 
stellar lines may be unseen in this low signal-to- 
noise region of the spectrum. 

The stellar atmosphere for this type of star also 
has broad Lya + He II absorption. It is, therefore, 
important to account for the stellar continuum un- 
derlying the interstellar absorption in our fitting 
process. For this purpose we used the synthetic 
spectrum and stellar Lya + He II line profile de- 
scribed in §3. 

The computed stellar spectrum was shifted to 
a heliocentric velocity of +20.2 km s -1 , the ra- 
dial velocity of BD +28° 4211, for comparison with 
the observed Lya profile. We determined this ra- 
dial velocity by measuring the velocities of promi- 
nent stellar lines in the STIS E140M observa- 
tions, the average of which gives a heliocentric 
velocity of Vhei — +20.2 ± 0.8 km s -1 . The 
quoted uncertainty is the 1 a standard deviation 
of the individual measurements about the mean 
and does not include the uncertainties in the STIS 
absolute wavelength scale, which probably corre- 
spond to ~ 1.6 to 3.2 km s -1 . This determina- 
tion agrees well with the stellar radial velocity 
of Vhei = +21.4 ± 1.7 km s -1 derived by Her- 
big (1999), where we use the mean and la stan- 
dard deviation of his individual measurements. 
For comparison, the average velocity of interstel- 
lar material along this sight line is v^ei ~ — 9 to 
-13 km s" 1 (stronger lines give more negative ve- 
locity components due to the presence of weaker 
components on the blue side of the interstellar 
profiles). The stellar radial velocity is well con- 
strained, although we found that offsets from this 
best value as large as 5 km s _1 have no effect on 
the derived H I column density. 
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The model stellar spectrum was scaled to best 
match the observed flux distribution using a 
second-order Legendre polynomial, the parame- 
ters of which are treated as free parameters in 
the fitting process. In general, the shape of the 
calculated stellar model far from line center fits 
the data relatively well. However, the adoption 
of a low-order polynomial correction to the syn- 
thetic stellar continuum allows for several possible 
systematic effects, including mis-matches between 
the observed spectrum and the model, slight re- 
distribution of flux in the scattered light removal 
process, stellar atmospheric variability, and other 
properties of the detector and observations that 
could change the absolute flux distribution on 
Angstrom scales. 

We also allowed for the effects of uncertainty in 
the flux zero point, treating it as a free parameter 
in the fit. This corrects for any small uncertainty 
in the background subtraction on scales of several 
Angstroms. The derived corrections to the flux 
zero point were always less than a percent or two 
of the continuum away from the line core. 

Because of the possibility that the stellar Lya 
profile could be variable, we fit separately each 
of the individual observations taken at different 
times. No evidence of variability was found. This 
is discussed in more detail below in our discussion 
of systematic uncertainties. 

4.4. Results and Discussion of Systematic 
Uncertainties 

Figure 4 shows the best fit intei stellar H I pro- 
file on top of the coadded STIS observations. The 
final estimate for the H I column is log TV (H I) = 
19.846 ^°q 3®, where the uncertainties are 95% con- 
fidence limits (2 a if Gaussian errors are appropri- 
ate) including systematic effects. This is an aver- 
age of the results from the two A (H I) determi- 
nations, which differed by only 0.003 dex (< 1%), 
and uncertainties, for which we adopted the larger 
of the uncertainty estimates. 

The uncertainties in A r (H 1) given above in- 
clude estimates of the statistical (random) and 
systematic uncertainties added in quadrature. Be- 
cause of the high quality of the data and the 
strength of the Lya damping wings, the system- 
atics dominate the statistical uncertainties in our 
analysis. Given the importance of the systematic 


uncertainties, we discuss several below and their 
impact on the A r (H I) determination. 

Stellar model uncertainties - The largest iden- 
tified source of systematic uncertainty for N(H I) 
is the adopted stellar model. However, fits to the 
observed Lya profile which use only a straight line 
continuum rather than a stellar model for nor- 
malization produce H I column densities that are 
larger by only +0.1 dex (~ 25%). This sets an up- 
per bound on the magnitude of the uncertainties 
due to the stellar models. 

The effects of imperfectly-known stellar atmo- 
spheric parameters on A r (H I) were estimated by 
refitting the interstellar Lya absorption profile us- 
ing model atmospheres with the most extreme 
stellar properties that are still nominally consis- 
tent with the determination of the atmospheric pa- 
rameters (see §3). For this purpose we chose those 
models with the strongest and weakest stellar Lya 
absorption profiles from a grid of atmospheres cov- 
ering a wide range in parameters: logHe/H(— 0.6 
to -1.4), T eff (78,000 to 86,000 K), and log g (6.0 
to 6.4). 

The N ( H I) values derived in this way were 
treated as extrema corresponding to 95% con- 
fidence limits (2a assuming Gaussian statistics) 
about the best-fit result. While this approach does 
not account for any uncertainties in our knowl- 
edge of the fundamental physics of stellar atmo- 
spheres or radiative transfer, the adopted extreme 
atmospheres encompass a sufficiently large range 
of parameter space that the effects of any sub- 
tle atmospheric physics are expected to be small 
in comparison to the uncertainties estimated in 
our approach. The uncertainties associated with 
the adopted stellar model and our allowance for 
stronger and weaker stellar Lya profiles dominate 
the uncertainties in the N(H I) determination. 

We note that the effects of line-blanketing 
within the model atmosphere will be much smaller 
than the effects of the extreme atmospheric pa- 
rameters investigated above. We verified this as- 
sumption with appropriate model calculations, 
and that the variation in the stellar Lya profile 
occurs primarily in the core of the line where the 
residual flux in the observed spectrum is zero. 

Residual scattered light - In principle, residual 
scattered light in the STIS E140M spectra could 
bias the N(R I) determination. The core of Lya 
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is indeed very near zero, and we allowed the true 
flux zero point to vary as a free parameter in the 
fits. However, if light were distributed in the line 
wings in an appropriate way during the scattered 
light removal process, this could potentially be a 
source of systematic uncertainty in our H I column 
density determination. 

As a check of the N ( H I) value derived 
from STIS E140M observations, we also analyzed 
archival observations taken with the pre-CO STAR 
Goddard High Resolution Spectrograph (GHRS) 
first-order G160M grating through the small sci- 
ence aperture (0725 x (y/25). These spectra have a 
resolution of ~ 15,500 (~ 19.3 km s” 1 FWHM). 
While the observations do not cover the entire ex- 
tent of the red wing of Lya, the holographically- 
ruled G160M grating had excellent scattered light 
properties (Cardelli, Ebbetts, & Savage 1993). 
Eight observations with a total exposure time of 
1152 s were reduced and combined using proce- 
dures described by Howk et al. (1999). From these 
spectra we found log A r (H I ) = 19.842 ± 0.020 (2a, 
statistical uncertainties only). This value is so 
close to the value derived using the STIS E140M 
observations that we believe residual scattered 
light uncertainties are not likely to be a signifi- 
cant source of error. 

Stellar atmospheric variation - The atmo- 
spheric profile of the star could be variable, and 
there is some evidence that BD +28° 4211 may 
be part of a binary system (see Napiwotzki 1999; 
Massey & Gronwall 1990). To account for the 
possibility that the stellar Lya profile might vary, 
we separately calculated the best-fit H I column 
densities for each of the eight individual STIS ob- 
servations. The dispersion in these best-fit mea- 
surements for the ensemble of observations was 
small and was added in quadrature to the final 
uncertainty estimate. The dispersion of the indi- 
vidual measurements about the mean was much 
smaller than that resulting from the variation of 
stellar model parameters described above. 

5. COLUMN DENSITIES OF ATOMIC 
DEUTERIUM, NITROGEN AND OXY- 
GEN 

5.1. General Considerations 

Both curve of growth (COG) and profile fitting 
techniques were used to determine N ( D I) and 


N ( 0 I) in order to better understand potential 
sources of systematic error. Profile fitting alone 
was used for N (N I) because there are only three 
unsaturated and unblended N I lines available (the 
A952 multiplet). The FUSE spectra were also an- 
alyzed to develop a model for interstellar N( H 2 ) 
for rotational levels J — 0 — 4. The compari- 
son of measurements from multiple channels (SiCl 
and SiC2) and observations (LWRS and MDRS) 
is a powerful tool that improves our statistics and 
helps to identify poor data quality due to instru- 
mental effects. 

In the COG method the equivalent width, W\, 
of several lines of the same species having different 
oscillator strengths are measured and their distri- 
bution in the \og(W\/\) - log(7V/A) plane is com- 
pared with a theoretical curve of growth for a sin- 
gle absorption component subject to a Maxwellian 
velocity distribution. The free parameters in the 
fit are the column density, N and Doppler spread 
parameter, b. However, if there is a complex ve- 
locity distribution, b is regarded as an “effective” 
Doppler spread and as such has no physical mean- 
ing. Continuum placement and equivalent width 
measurements are treated in independent steps 
and assume some a posteriori conservative error 
bars to take into account the uncertainties in con- 
tinuum placement. The COG method is poten- 
tially vulnerable to error because of the assump- 
tion of a single- component profile. In particular, 
at the FUSE spectral resolution narrow, unre- 
solved components could be concealed by broader 
features, which could result in underestimated col- 
umn densities (Nachman & Hobbs 1973, Savage & 
Sembach 1991). Jenkins (1986) showed that col- 
umn densities accurate to ~10% can be obtained 
with the COG method provided that the lines are 
not heavily saturated (tq < 2) and have smooth 
distributions of b and r in their components. In 
the case where all transitions are on the linear part 
of the COG, the column density is insensitive to 
the line-of-sight velocity structure. One of the ad- 
vantages of the COG method is that it is not too 
sensitive to the LSF. On the other hand, transi- 
tions blended with lines from other species should 
be avoided in this analysis technique. In this case, 
profile fitting may be more appropriate. 

The profile fitting method makes the same basic 
assumptions as the COG: an interstellar absorp- 
tion line is formed in a single component and is 
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described by two parameters: N and b. In ad- 
dition the wavelength centroid of the line, or ra- 
dial velocity, v ry is also a free parameter. The 
fitting technique also assumes thai. each interstel- 
lar component produces absorption lines described 
by Voigt profiles. We note that in our analysis the 
Doppler parameter, 6, has no unique physical in- 
terpretation. Normally b 2 = 2kT/m-\-v^ urb . How- 
ever, from the FUSE spectra alone, with spectral 
resolution of ~ 15 km s -1 , tempe rature (T) and 
turbulent motion (u* ur .b) effects on the line widths 
are not easily separated. We can, however, use b 
to set an upper limit on T. 

Lines of sight with multiple c omponents are 
often more reliably analyzed with a fitting tech- 
nique than it is with the COG. In practice, to fit a 
multiple-component COG requires a prior knowl- 
edge of v T and b for each component as well as 
their expected relative ratios. This information is 
not necessary with profile fitting be cause in princi- 
ple it lies in the profiles. For the same reason, the 
blending of adjacent lines is also easier to model 
with a fitting procedure. 

A second significant advantage of our pro- 
file fitting analysis is that absorption lines and 
continua are fit simultaneously so that contin- 
uum shape and placement are treated as free pa- 
rameters. We have used the profile fitting pro- 
gram Owens. f (Lemoine et al. 2001, Hebrard 
et al. 2001). Owens. f is a powerful profile fit- 
ting program that simultaneously models multi- 
ple spectral lines, species, and absorption com- 
ponents in an arbitrary number of spectral re- 
gions, or windows. Many quantities characteriz- 
ing the fit may be treated as free or fixed param- 
eters (N,T, Vturb, LSF, continuum placement and 
shape, background level, velocity offset between 
spectral windows, etc.). The program uses an op- 
timized x 2 ‘ m i n i m i za ti° n algorithm to obtain fast 
convergence. The flexibility of Owens . f allowed 
us to analyze many different cases. As with any 
convergence algorithm, parameter space was ex- 
plored to verify that a converged solution was the 
true minimum and not a local minimum far from 
the true solution. This was done by changing the 
initial parameters and/or conditions and then run- 
ning the program again to reach a new converged 
solution. A potential disadvantage of profile fit- 
ting is that knowledge of the instrument LSF is 
required, and errors in the assumed LSF have the 


potential to bias the results. 

The instrumental LSF is a key parameter in 
the profile fitting analysis of interstellar absorp- 
tion lines. The characterization of the FUSE LSF 
was incomplete at the time of this analysis (see 
Kruk et al. 2001). Although there is evidence 
that the FUSE LSF has wings that effect the zero 
level of heavily saturated lines (§2), two indepen- 
dent analyses (Hebrard et al. 2001, Wood et al. 
2001) show that modelling of weak unsaturated 
lines with a single Gaussian LSF and a double 
Gaussian LSF give very consistent results within 
the uncertainties of the FUSE LSF parameters 
(see the companion papers referenced above for 
further details). For our profile fitting analyses 
we used a single gaussian LSF for each component 
modelled. 

In the next subsection the interstellar H 2 model 
is derived. The following subsections describe the 
COG analysis of D I and O I and the profile fitting 
analysis of D I , O I, and N I. 

5.2. Interstellar H 2 Column Densities 

The FUSE spectra of BD +28° 4211 contain a 
large number of absorption lines from interstellar 
molecular hydrogen. Absorption arising from H 2 
rotational levels J — 0 — 4 can be clearly seen 
along this line of sight. Many H 2 lines axe blended 
with stellar and other ISM features. For transi- 
tions corresponding to the rotational level J = 4 
there are not enough unblended lines to derive a 
reliable column density using either COG or pro- 
file fitting techniques. However, a 3cr upper limit 
for N(H 2 ) J — 4 is derived. A few H 2 transi- 
tions from J = 5 and 6 may also be present. The 
best candidate is Werner 0 — 0 Q(5) at 1017.831 
A. However, the FUSE spectra of BD +28° 4211 
contain many unidentified stellar features and it 
is possible that this and other weak transitions 
could be stellar in origin, or if actual H 2 ISM lines, 
blended with stellar lines. Due to these uncertain- 
ties, the column densities for J > 5 could not be 
estimated. The H 2 column densities for rotational 
levels J = 0 — 3 were measured using both pro- 
file fitting and COG methods, using lines that ap- 
peared free of blending, as judged by the shape 
of the absorption profile and with the help of the 
stellar model. 

Profile fitting used only unsaturated lines, i.e. 
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lines for which the absorption profile deconvolved 
with the LSF did not have a residual intensity 
in the line core < 0.1. The column density 
for each rotational level was determined indepen- 
dently. The synthetic stellar spectrum gave guid- 
ance when selecting the H2 lines to be analyzed 
in order to avoid blending with stellar features. 
Each spectral region (several Angstroms) was first 
normalized with a spline function. A subset of 
this region, a spectral window of ~ lA, was used 
by Owens . f where a polynomial (up to i th order) 
was fit to the continuum away from the line. We 
restricted the size of the spectral windows used 
in the Owens . f analysis because of concern about 
small nonlinearities in the FUSE wavelength cali- 
bration on larger scales. The number of lines an- 
alyzed for each J level are listed in Table 3. All 
the H2 lines analyzed lie away from the H I lines. 
The background was set to zero for each window. 
The width of the Gaussian LSF was free to vary 
from window to window and between successive 
fits. Once the best fit was found, corresponding 
to a minimum x 2 , the error in the column density 
was calculated by fixing N(R 2 ) to a different value 
and computing a new optimized fit, which always 
had a larger x 2 - The difference of the new reduced 
X 2 from Xmin determined the 1, 2, 3, 4, 5, 6, and 
7a deviations from the best fit. The adopted la 
error was the average of these values to allow for 
any asymmetry in the Ax 2 curve. Only an upper 
limit was derived for N( H2) J = 4, the result of 
its weak transitions being indistinguishable from 
weak stellar features. jV(H 2 ) and 2 a errors de- 
rived this way are given in Table 3. 

The H2 equivalent widths were measured by di- 
rectly integrating the area of several lines for each 
rotational level for J — 0 — 3. The continuum near 
each line was defined by a Legendre polynomial up 
to 5 th order. W\ of each transition is the weighted 
average of W\ measured in each channel where 
the line was present. A single-component Gaus- 
sian COG was fit to the data for each rotational 
level, and the residuals about the best-fit curve 
were minimized (Figure 5). The resulting column 
densities are completely consistent with the profile 
fitting results, although with somewhat larger un- 
certainties. Assuming the population of the lower 
rotational states of H 2 is determined by collision al 
excitation and a Boltzmann distribution, the exci- 
tation temperatures are Ti )0 = 313 K , T 2( i — 294 


K and T 3,2 — 296 K. However, the similarity of T{ t j 
for J = 0 — 3 indicates that one or more of these 
assumptions may not be correct. The uncertainty 
in N(R 2 ) J = 0 may permit a more typical value 
of Ti 0 ~ 100 K. The total H 2 column density is 
log(7V(H 2 ))= 15.13^o os cm" 2 . 

5.3. D I and O I Curve-of- Growth Analy- 
ses 

The high S/N of the FUSE observations of 
BD 4-28° 4211 permits analysis of the D I tran- 
sitions down to 916 A where the H I profiles start 
to overlap. However, some of the higher D I Ly- 
man series lines (916-918 A) have strengths that 
vary significantly from line to line and channel to 
channel. These inconsistencies probably originate 
in the fixed pattern noise of the detector and/or 
in blends with weak stellar features or interstellar 
H 2 . Therefore, the COG analysis of D I was lim- 
ited to the six transitions listed in Table 4. Most 
of the D I lines (Ly7 - Ly 77) lie in the blue wing 
of the corresponding H I lines. Four out of the six 
transitions fall on the linear part of the curve of 
growth. Hence, the velocity structure on the line 
of sight does not influence the final column density 
estimate (see §5.4). 

5.3.1. Dl Continuum Normalization 

The stellar H I and He II and interstellar H I 
absorption profiles dominate the shape of the con- 
tinuum near the D I feature. A continuum model 
was constructed, where the key parameters are 
N(R I) , JV(H 2 ) , and b for the ISM model com- 
puted with Owens. f and T e //, logg and He/H for 
the stellar model. We adjusted the stellar plus in- 
terstellar model to fit the data (a small zero-point 
shift in the model’s wavelength scale and scaling 
the absolute flux of the model to match the ob- 
served flux level - see Fig. 3). These are discussed 
further below. 

The stellar model has been discussed in detail in 
§3. The uncertainties in T e // and log g considered 
in the N ( H I) analysis (~5%) have a negligible 
effect on the shape of the HI + He II profile near 
the D I lines we are studying. While varying these 
parameters would change the goodness of the fit 
far from the core of the D I lines, the continuum 
placement would be similar. As to the He/H ratio, 
this affects only the two strongest lines (D I Ly7 
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and Ly£) where He II is readily apparent in the 
model. 

The overall shape of the model agrees well with 
the spectrum, although we noted some localized 
discrepancies. The different spect ra showed flux 
scale differences from the model of up to +5% for 
LWRS and +10% for MDRS. The latter suffered 
from small variations in slit transmission due to 
mirror motion (see Sahnow et ah 2000) and a 
small part of the flux was lost in this process. The 
MDRS flux calibration is also less well known than 
that of the LWRS channel. 

We modelled the interstellar H l profiles by us- 
ing Owens. f to fit the relevant H I Lyman se- 
ries transitions (919 - 972 A), for< ing the column 
density to log(7V(H I) ) = 19.85 and the adopted 
the 7V(H 2 ) parameters from Table 3. A one- 
component fit agreed well with the data. A three- 
component fit was also performed (cf. §5.4) but 
there was no improvement in the fit because in 
H I these components are buried in the saturated 
line cores. The purpose of this H 1 model was to 
find a good empirical match with the data to assist 
in the measurements of D I W \ , n :>t to determine 
N ( H I) . 

5.3.2. W\ Measurements 

The general method for measuring W\ for the 
D I lines is illustrated in Figure 6. First, the data 
were divided by the stellar model appropriately 
registered in A and F\ (Fig. 6a and d). This 
process sometimes produced what was judged to 
be a poor fit in the red wing of interstellar H I, 
especially for D I A925. We believe this is due to 
differences in the number, location, and strength of 
weak stellar lines between the blue and red wings. 
In spite of the concern about such line blanketing 
differences, W\ for D I A925 is consistent with the 
other D I transitions analyzed. 

Second, the normalized spectrum was then cor- 
rected for interstellar H I (Fig. 6b and e) by us- 
ing the H I model profiles described above. W\ of 
each D I feature was measured by fitting a Gaus- 
sian to the residual D I absorption line (Fig. 6c 
and f) using ail of the profile except the far red 
wing of where the near-zero flux n<*ar the line core 
of H I produces large residuals. Since the prin- 
cipal source of error in this analysis is the stellar 
continuum placement, we adopted a conservative 


approach, scaling the continua by ±3% and re- 
measuring W\. This produced noticeably poor fits 
that represent extreme cases. We assumed that 
these fits were 2cr from the best fit. Therefore, the 
JV(D I) errors include the systematic error due to 
continuum placement. 

Table 4 gives the D I W\ values for each chan- 
nel and aperture. There are four independent 
measurements of W for each DI transition (SiCl 
and SiC2, LWRS and MDRS). These measure- 
ments were consistent to within la, except for the 
three cases where W\ was omitted. These mea- 
surements were discarded because they were in- 
consistent (discrepant by >50%) with the other 
measurements of the same transition. We con- 
cluded that these cases were probably caused by 
detector artifacts (fixed-pattern noise) after care- 
ful comparison of line profiles in the four aperture 
and channel combinations. The far right hand 
column in Table 4 is the mean W\ used in the 
COG. The COG fit for D I is shown in Figure 
7. We found log N{ D I) = 14.99 ± 0.10 (2a) and 
b ~ 5.2 km s -1 . 

5.3.3 . O I Curve of Growth 

A curve of growth was constructed for 0 I, 
for each segment (SiCl and SiC2, MDRS and 
LWRS), using the lines and /-values of Table 4. 
O I AA924.52, 936.63, and 950.88 were not used 
because of blending with stellar or H 2 features. 
The general procedure to construct the curve of 
growth is similar to what was described above for 
the H 2 COG analysis (§5.2). For the O I COG 
we chose lines that are free of blending with other 
ISM or stellar features. Unlike D I and H 2 , O I /- 
values may have relatively higher uncertainties in 
the FUSE bandpass, all of them coming from the- 
oretical calculations (Morton 1991; Morton 2000, 
private communication; Verner, Barthel, & Tytler 
1994). Unlike the pathological case of O I A1026.5 
pointed out by Jenkins et al. (2000), the other 
O I /-values appear to be self-consistent. W\ 
and associated errors measured for the four spec- 
tra are listed in the Table 4. The COG fit to the 
mean W\ is shown in Figure 8 and the derived 
column densities are given in Table 5. A COG 
fit to the weighted average of W\ for the different 
channels and apertures yields log N ( O I) cog — 
16.23 ± 0.08 cm -2 (2a) and b = 4.5 km s" 1 . 


11 



5.4. Profile Fitting Analysis 

A profile fitting analysis of D I , N I, and 0 I 
column densities wets performed to provide a com- 
parison with the results of the COG analysis in 
the previous subsection. The general technique is 
identical to that discussed in §5.2. Since Owens. f 
was used to fit D I , N I, and 0 I simultaneously, 
the following discussion also applies to N ( O I) and 
iV(N I). 

The five D I transitions analyzed (Lye, Ly£, 
Ly 77 , Ly0, and Ly l) are not saturated. With the 
use of a fitting procedure, the choice of the lines to 
analyze is different than with the COG technique. 
An example is Ly£ which was not used in the COG 
analysis because it is blended with H 2 J = 1 and 
J — 2. By adopting the H 2 solution (Table 3) as 
a constraint in the fit, this transition is now part 
of the fitting analysis and easily deblended from 
the H 2 . We note that the errors in JV(H 2 ) are 
not significant for the N ( D I) analysis because 
the H 2 lines blended with Ly£ are on the flat part 
of the COG. The two strongest D I transitions 
used for the COG (Ly 7 and Ly 5) were omitted 
from the profile fits because they are mildly satu- 
rated (t 0 ~ 2). Uncertainties in our knowledge of 
the FUSE LSF could influence the results derived 
from profile fits to such lines. In any case, N(D I) 
is well constrained by the other five transitions. 

Although the continuum close to the O I lines is 
very smooth and easily fit with a low-order polyno- 
mial, some of the O I transitions used in the COG 
are indeed saturated and for this reason they were 
not included in the N ( 0 I) profile fitting analysis 
using Owens. f. Only three O I transitions were 
used to determine N ( O I): 919.92 A, 925.45 A, 
and 930.26 A. The strongest of these has r 0 ^ 1. 

There are many N I transitions in the FUSE 
range and four multiplets are detected : AA952, 
AA953, AA964, and AA1134. The 952 A triplet 
is free of any blending with H 2 or ISM atomic 
lines and seems clear of stellar features. Since 
all the lines of this triplet are still on the linear 
part of the curve of growth, deriving the column 
density is straightforward. This is not the case 
for the AA953, AA964, or AA1134. These multi- 
plets ets are all saturated and much less sensitive 
to N than they are to 6. The 964 A triplet is 
the next strongest after AA952, but these lines are 
blended with several H 2 lines, P II A963.8, and 


with unidentified stellar features. With only one 
multiplet (AA952) we did not perform a COG anal- 
ysis of N I. The N I /-values used were provided 
by D. Morton (2000, private communication). 

The Owens . f analysis of the D I , N I, and 0 I 
lines allow many quantities that characterize the 
fits to be treated as free parameters. These include 
the wavelength and flux zero points of each profile, 
the polynomial coefficients for each spectral win- 
dow, the LSF parameters (discussed below), the 
column densities N , and Doppler parameters b . In 
total, there were about 40 spectral windows and 
about 1500 degrees of freedom to fit 11 transitions. 
The H 2 column densities (Table 3) were adopted 
as a fixed constraint. We fit the five lines of D I , 
three of O I, and three of N I in the four data sets 
simultaneously, providing a global solution. 

The continuum for the D I line profiles was de- 
termined by a different approach than was used 
for the D I W\ measurements. There is sufficient 
information immediately adjacent to the D I fea- 
tures to establish the continuum with a polyno- 
mial without normalizing the spectrum or mod- 
elling the stellar and interstellar H I . The shoul- 
der of the blue wing of the interstellar H I profile, 
excluding the D I feature, was fit by a i th order 
polynomial. Examples of the continuum and D I 
profile fits are shown in Figure 9. The effective b 
is not tightly constrained in this analysis because 
the true shapes of the D I profiles are unresolved 
for these weak lines. 

The effect of simultaneously modelling the in- 
terstellar H I and D I profiles and fitting a 
polynomial to the stellar contribution alone was 
also examined. The line-of-sight velocity struc- 
ture becomes critical when higher Lyman series 
of D I and H I are fit together. In this case for 
BD +28° 4211 the effective b is dominated by H I 
and is larger than the D I effective 6, much larger 
than the expected y/2 difference between H and 
D. This is probably the result of very weak H I 
components not visible in D I or any other species 
that become a factor in the width high Lyman se- 
ries profiles of H I . Overestimating b for D I leads 
to an anomalously low N(D I) . When D I and H I 
are fit together a larger fitting window (AA ~4 
A) is also required, forcing D I to be shifted from 
its true center as a result of small residual nonlin- 
earities in the FUSE wavelength calibration. Due 
to the introduction of these systematic errors we 
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rejected this fitting approach in favor of the tech- 
niques described in the previous paragraph and 
illustrated in Figure 9. 

The influence of the LSF assumptions and the 
velocity structure were carefully examined. Three 
cases were considered in the proiile fitting anal- 
ysis. One way of investigating the sensitivity of 
the fit to the assumed LSF was to let Owens./ ad- 
just it as part of the x 2 minimization. For Case 
1 we considered a fixed LSF associated with each 
channel /aperture combination, ignoring the possi- 
ble variations with wavelength. The choice of the 
LSFs was made after several Owens . f test runs 
with different parameters. In Case 2, we assumed 
a LSF modeled by a single Gaussian with FWHM 
that varied freely with wavelength and from chan- 
nel to channel. Case 3 examined the effect of line- 
of-sight velocity structure as derived from STIS 
Fe II profiles on the derived column densities. The 
continuum was determined in all three cases by 
the technique illustrated in Figure 9. For all three 
cases, the LSF was modelled as a single Gaussian 
for reasons discussed at the end of §5.1. 

A free LSF (Case 2) clearly improved the fit, as 
indicated in Table 6, compared with the fixed LSF 
(Case 1). Only a few spectral windows showed 
LSF width variations as large as 20% (< 2 pixels). 
Such variations may be the result of fixed pattern 
noise that is only partly corrected by means of 
the FP-split procedure, or from weak, unresolved 
photospheric features. For Case 2 the reduced \ 2 
was slightly lower than before but the 7V(D I) and 
N ( O I) were unchanged. The fit of the N I lines, 
fit simultaneously with D I and 0 1, showed slight 
differences but are within the 2cr uncertainty. Al- 
though the LSF calibration as a function of A, 
channel, and aperture is not yet well determined 
(Kruk et al. 2001), we accounted for this by con- 
sidering two extreme assumptions (free vs. fixed 
LSF) and adopting an average column density that 
includes this systematic effect in the errors. 

For Case 3 we used high resolution STIS E230H 
spectra to provide resolved Fe II AA2344, 2382 
line profiles, the only Fe II lines available in 
the existing archival high resolution spectra of 
BD +28° 4211. Profile fits to the STIS spec- 
tra alone found three components at Vhei — 
— 2.4, —9.7, and —17.0 km s _1 (listed in order 
of decreasing Fe II column density) and a total 
Fe II column density of log 7V(Fe II)= 13.81 cm -2 


(see Figure 10). We also performed a simultane- 
ous fit of FUSE and STIS data using Owens. f. 
This fit found the velocity separations of the 
three components to be 3.5 and 10.9 km s -1 
and log A r (Fe II)— 13.79 cm -2 , compared with 7.3 
and 7.3 km s -1 for the STIS data alone. These 
differences (3 — 4 km s _1 ) reflect the level of un- 
certainty of these fits. 

For the Case 3 fits, the relative ratios of each 
component for the neutral species was very dif- 
ferent from that in Fe II. Specifically, the weakest 
(bluest) Fe II component had no contribution from 
D I , N I, or O I. This could be due to ionization 
effects or iron might be depleted in different ways 
in different environments, while N I and O I are 
only lightly depleted (Savage & Sembach 1996). It 
appears that the fit of the neutral species suffered 
from the use of the Fe II lines, showing a higher re- 
duced x 2 - We believe that the neutral species are 
not fit well by the Fe II velocity structure because 
the 6- values and/or the velocity structure appears 
to be different. Since we do not resolve the com- 
ponent structure with FUSE , we rely on the one 
component solution for D I , N I, and O I. 

The column density errors were calculated by 
computing new fits for different fixed values of N 
away from the best fit column density and deriving 
Ax 2 as a function of N. The 3a error is given by 
Ax 2 = 9. The average gives log N(D 1) = 14.99 i 
0.03 cm _2 (2a). 

The COG and profile fitting analyses give the 
same D I column density and have error estimates 
that differ only slightly at the 2 a level. It is reas- 
suring that both approaches give results within 5% 
of each other. The errors from the COG analysis 
are more conservative than the ones coming from 
the profile fits. The reason for this is that mov- 
ing the continuum by ±3%, as was done for the 
COG, we have probably overestimated the errors. 
Yet they seem to fairly reflect the scatter seen in 
the individual W\ measurements. Combining the 
results, we adopted a straight average of N( D I) 
from the COG and profile fitting analyses and ar- 
rive at the final D I column density log N( D I) 
= 14.99 ± 0.05 cm -2 (2cr) for the BD +28° 4211 
line of sight. 

Figure 11 compares the O I fits for SiC2 spec- 
tra from the LWRS and MDRS observations. 
Combining the results from the different LSF 
approaches (see Table 6) gives log N(0 I) = 
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16.21 ± 0.10 cm 2 (2cr ) from the profile fitting 
analysis. 

Most of the 0 I lines used in the profile fitting 
analysis are close to saturation which may make 
both this and the COG method sensitive to the 
velocity structure. We investigated this possibil- 
ity by analyzing the effect of a three-component 
velocity structure on the COG analysis. The goal 
of this investigation was to see if a one-component 
COG could still be applied with success on this 
line of sight. We computed a model for the line 
of sight for D I and 0 I lines with three compo- 
nents as seen in the STIS Fe II profiles. Perform- 
ing a one-component COG analysis of the model 
line profiles yielded D I and O I column densi- 
ties within 5% of the model and an average ef- 
fective b of 5.8 km s'" 1 . Although the velocity 
structure found in Fe II is not identical to that for 
the neutral species, it is the best available. For 
the BD ±28° 4211 line of sight a one-component 
COG gives reliable results. We find that the ef- 
fective b from the analysis of the three-component 
simulated profiles is consistent with the effective 
b from the analysis of the actual data. The final 
value adopted for the O I column density is log 
N(0 I) = 16.22 ± 0.10 cm" 2 (2a). 

Examples of the fits to the N I AA952 lines are 
shown in Figure 12 and the results listed in Ta- 
ble 6. The mean fit gives log N(N I)= 15.55 ± 
0.13 (2a). 

6. RESULTS AND DISCUSSION 

The abundances derived in the previous sec- 
tions are summarized in Table 8 along with their 
ratios. The 2 a errors in the abundance ratios 
are those from the individual column densities 
combined in quadrature. We find that D/H = 
(1.39 ±0.21) x 10~ 5 (2a) toward BD +28° 4211. 
Golz et al. (1998) observed BD +28° 4211 with the 
Tubingen echelle spectrograph (A/ A A ~ 10,000) 
on the ORFEUS-SPAS II mission and found D/H 
— (0.8^0 4) x 10 _5 (la). Our result is at the up- 
per end of their range. The FUSE result for 
BD ±28° 4211 is just marginally consistent with 
the higher D/H value found for the more distant 
sight line toward Feige 110 (D/H — (2.1 ± 0.70) x 
10 -5 ) by Friedman et al. al. (2001). The to- 
tal neutral gas column toward BD +28° 4211 is 
very similar to that toward j 2 Vel, C Pup, an d 


S Ori A (at distances of 250-500 pc), the sightlines 
studied with IMAPS to measure D/H (Jenkins et 
al. 1999, Sonneborn et al. 2000). D/H toward 
BD +28° 4211 is in the middle of the range of val- 
ues found in the IMAPS studies. 

Our result for D/H toward BD +28° 4211 is 
similar to values found on lines of sight in the Lo- 
cal Bubble by FUSE (Moos et al. 2001) and HST 
(Linsky et al. 1998), even though BD-f28°4211 
lies just outside the Local Bubble. The total 
H I column density toward BD ±28° 4211 is log 
jV(H I) = 19.85, much higher than the other LISM 
sight lines studied thus fax by FUSE for D/H (ex- 
cept for Feige 110). The Na I contours for the 
LISM (Sfeir et al. 1999) show that BD +28° 4211 
is located in a region of higher column density just 
beyond the edge of the Local Bubble, with much 
lower densities in the intervening ~ 100 pc. The 
H I column toward BD +28° 4211 is large enough 
that the local cloud is a negligible contribution. 
Our abundance measurements are therefore dom- 
inated by the gas ~ 100 pc from the Sun. As 
discussed by Moos et al. (2001), there is now 
considerable evidence for a constant D/H within 
the Local Bubble. Beyond ~ 100 pc, however, or 
more specifically outside the Local Bubble and log 
N ( H I) ge 19.5, there appears to be a considerable 
spread in D/H as larger amounts of gas are probed. 

The total H2 column density is ~ 1.3 x 

10 15 cm -2 , and the molecular fraction /(H2) ~ 
3.8 x 10” 5 . This low value of /(H 2 ) is typical of 
lower H I column density ISM sight lines (Tum- 
linson et al. 2001). The abundance of H 2 at 
~ 100 pc indicates that some molecular gas is 
present, but the low column densities and high 
excitation temperature are consistent with a low 
density environment where UV pumping is prob- 
ably significant. 

The profile fitting and COG analyses for D I 
and O I consistently found an effective 5-value of 
~ 5 km s _1 . Although FUSE lacks the spectral 
resolution to sepatate thermal and non-thermal 
effects on the neutral gas line profiles, we can 
set an upper limit to the temperature. For D I 
b = 5.0 km s" 1 corresponds to T < 6000 K. For 
O I, however, the same b corresponds to T <40,000 
K. There must either be a significant non-thermal 
contribution to the line broadening of the neu- 
tral species, or the component structure on the 
BD +28° 4211 line of sight has a significant effect 
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of the line shapes. Very high resolution spectra of 
species such as Na I are needed to determine the 
component structure, T and Vt. ur , of the neutral 
gas as a function of velocity. 

The O/H ratio toward BD +28' 4211 is (2.37 ± 
0.55) x 10“ 4 (2 a), a value that is anomalously low 
with respect to the Local Buble and the ISM. The 
BD +28° 4211 O/H is marginally consistent with 
the mean O/H ratio in the ISM (<0 /H>/sa/) 
of (3.19 ± 0.28) x 10~ 4 (2a) found by Meyer et 
al. (1998) in a study of O I A1355 on 13 sight 
lines toward OB stars at distances of 130 - 1500 
pc. Recently, however, Meyer (2001) updated 
<0/H> jsm to (3.43 ± 0.28) x 1 0 ~ 4 as a result on 
a revised /-value for the A 1355 transition. The 
FUSE O/H ratio toward BD 4-28° 4211 is well 
outside the estimated uncertainties in the revised 
<0/H>jsm- In the Local Bubble <0 /H >lb= 
(3.9±0.3) x 10” 4 (la) (Moos et al 2001), slightly 
higher than Meyer’s revised value for the ISM, al- 
though consistent within their respective uncer- 
tainties. Moos et al. (2001) also finds the same 
value of the mean O/H when the more distant 
sight line to Feige 110 is included. Of the five 
sight lines for which O/H is measured by FUSE , 
only BD +28° 4211 gives a value that disagrees 
with <0/H> ism and <0 /H>x /j b If correct, the 
gas just outside the Local Bubble has a low O 
abundance. This implies that the gas toward 
BD +28° 4211 beyond the Local Bubble would ap- 
pear to have a different chemical history than that 
sampled in the Meyer et al. study 

We have searched for other possible explana- 
tions for the low O/H ration toward BD +28° 4211. 
O I is not vulnerable to ionization effects (Jenk- 
ins et al. 2000), so the O I-to-H I ratio should 
accurately reflect the total O/H abundance ra- 
tio. If the error is in the H I column den- 
sity, then correcting O/H to obtain agreement 
with <0 /H>/sm would result in a large D/H 
ratio (~ 2.0 x 10“ 4 ). Such a D/H value for 
BD +28° 4211 would be similar to that found for 
Feige 110 (Friedman et al. (2001) and q 2 Vel (Son- 
neborn et al. 2000). However, a ~ 40% increase in 
iV(H I) would be inconsistent wh h the available 
spectra. 

The D/O ratio for BD +28° 4211 appears to be 
high (Moos et al. 2001, Hebrard et al. 2001), argu- 
ing for N(0 I) as the anomalous quantity. Agree- 
ment with <0/H >jsm would require a ~ 40% 


increase in A r (0 I) . One or more of the O I lines 
analyzed might actually be saturated, even though 
our analysis shows (Fig. 11) that the three O I 
lines are unsaturated. Line saturation, possibly 
as the result of narrow (cold gas) components 
would result in an underestimated column density. 
The effect of potential line saturation in N ( O I) 
was examined by excluding the strongest O I line 
(A930) from the fitting analysis, leaving only the 
two w r eaker O I lines. For this case we found log 
N(0 I) = 16.28, which is within 2a of the best fit. 
We cannot determine whether this difference in N 
is the result of saturation, fixed pattern noise, LSF 
uncertainties, component structure, or some other 
cause. If this were the correct O I column den- 
sity, the O/H ratio would be 2.72 x 10 -4 , still well 
below < 0 /H>/ 5 m and <0/H >lb- 

The N/H ratio toward BD +28° 4211 is (5.08 ± 
1.66) x 10~ 5 (2a), a value that is marginally con- 
sistent with the mean ISM value of N/H = (7.5 ± 
0.8) x 10 -5 (2a) found for six OB stars by Meyer 
et al. (1997). However, the N II column den- 
sity could be significant (see below). If so, then 
the N/H ratio toward BD 4-28° 4211 could well be 
consistent with the mean ISM result. 

In the local ISM, the abundance of N I has been 
found to be reduced by ionization effects on some 
sight lines (Jenkins et al. 2000, Moos et al. 2001). 
N I has a photoionization cross section larger than 
that of H I and a large fraction of it can be ionized. 
There are only two N II transitions in the FUSE 
bandpass (A915.6 and A1085.5) and both of them 
are heavily saturated on all but the lowest column 
density sight lines. We derived a lower limit of log 
N ( N II) > 14.2 using both optical depth technique 
and profile fitting. Since we know that some N II 
is associated with the neutral component, we esti- 
mate an upper limit assuming that the b- value for 
N II is greater than that of the neutral gas. These 
assumptions give an upper limit of log 7V(N II) 
< 15.7. We cannot rule out the possibility that 
an apreciable fraction of the nitrogen on the sight 
line toward BD 4-28° 4211 could be singly ionized. 

Objects like the hot subdwarf BD 4-28° 4211 are 
important targets for studying D I and O I in the 
ISM with FUSE because they can sample regions 
of space beyond that accessible with white dwarfs 
( d < 100 pc) and closer that that sampled by the 
lightly-reddened O stars observable by FUSE ( d > 
1 kpc). Most of the nearer O stars far exceed the 
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FUSE brightness limit. 

High-quality spectra of the type shown in Fig- 
ures 1 and 2 are needed to understand these ob- 
jects, and the observational techniques to obtain 
them with FUSE are now available. Over time 
issues like stellar continuum placement and stel- 
lar line identifications may be better understood, 
allowing even more precise interstellar abundance 
measurements. 

This work is based on data obtained for the 
Guaranteed Time Team by the NASA-CNES-CSA 
FUSE mission operated by the Johns Hopkins Uni- 
versity. Financial support to U. S. participants has 
been provided in part by NASA contract NAS5- 
32985 to Johns Hopkins University. Support for 
French participation in this study has been pro- 
vided by CNES. 
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Velocity (km/sec) 


Fig. 2. — The higher Lyman lines in BD +28° 4211 from the SiCl LWRS spectra are shown on a velocity 
scale. The zero point is set to the observed center of the D I feature. The principal interstellar atomic 
transitions are indicated above th< spectrum and stellar transitions below. 
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Fig. 3. — Comparison between a portion of the FUSE spectrum (solid line) and a synthetic spectrum (dashed 
line) computed with T e ff = 82,000 K, log(He/H) = —1.0, logg = 6.2, and incorporating solar nitrogen and 
oxygen abundances. The photospheric lines are labelled at the top of the figure. The ISM lines are labelled 
just above the spectrum. 
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Fig. 4. — The wavelength region surrounding interstellar Lya towards BD +28° 4211. The histogram shows 
the coadded STIS E140M observations of BD +28° 4211. The dashed gray line shows the adopted stellar 
model ( T e ff = 82, 000 K, log# = 6 2, and logHe/H = —1.0) modified by a second order Legendre polynomial 
to bring it into agreement with the observed flux distribution. The two narrow absorption components in 
the stellar profile are the cores of FI I and He II. The solid line shows the best-fit H I model times this stellar 
continuum. The best-fit interstellar parameters are Vhei — — 13.3±0.8 km s _1 and log 7V(H I) = 19.846io“o36 
(both 2<t uncertainties). 
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Fig. 5. — Curve of growth for H 2 with fits corresponding to b = 2.8 km s -1 (dashed line) and b = 8.5 km s _1 
(dot-dashed line), which were the minimum and maximum b values found in the fit of the individual rotational 
levels. The la error bars are comparable to the scatter in the W\ measurements. They were omitted from 
the figure for clarity. 
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Velocity (km s _1 ) 


Fig. 6. — Continuum normalization for D l W\ measurements. Panels a) - c) illustrate the method used 
to measure W\ of D I A925.974 (Ly?y) in the SiC2 channel and LWRS aperture. Panel a) shows the stellar 
model placement (smooth, solid line) over the data (histogram-style line). Panel b) shows the ISM model (H I 
-P H 2 only) relative to the spectrum normalized by the stellar model. Panel c) shows the ISM-normalized 
spectrum from b) and the Gaussian fit to the D I Lyry profile used to measure the area of the line under 
the continuum. The red side of the reconstructed D I profile is affected by large residuals from dividing by 
the near-zero flux in the core of II I. Panels d)-f) show the same method applied to D I Ly 7 in the same 
channel and aperture. The weak f- ature in the blue wing of the stellar Ly 7 model profile (panel d) is He II, 
while the stronger feature directly below it in the observed spectrum is interstellar D I . 
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Fig. 7. — The best fit curve of growth for D I (solid line) corresponding to log N(D I) = 14.99 ± 0.10 (2a) 
and b = 5.2 km s” 1 . The dashed lines are the la errors in the fit. The equivalent width data points, shown 
with the wavelength of the transition and la error bars, are the weighted average of W\ from the SiCl and 
SiC2 channels and the LWRS and MDRS observations. 




log Xf 


Fig. 8. — The best fit curve of growth for 0 I (solid line) corresponding to log iV(0 I) = 16.23 ± 0.08(2a) 
and b — 4.46 km s” 1 . The dashed lines are the la errors in the fit. The equivalent width data points, shown 
with the wavelength of the transition and la error bars, are the weighted average of W\ from the SiCl and 
SiC2 channels and the LWRS and MDRS apertures. 
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Fig. 9.— The fits of D I Lye A937.5, LyC A930.5, Ly rj A926.0, and Lyt A920.7 for the SiC2 channel and LWRS 
aperture are shown on a relative velocity scale. The histogram lines represent the data, the thick solid lines 
are the final fits to the spectrum, and the dashed lines are the fits after deconvolution with the LSF. Note 
that no D I lines stronger than Lye were used since this transition is already optically thick (to ~ 2). Note 
also D I LyC is blended with H 2 (J = 1) and H 2 (J = 2). The profile fitting technique recovers the information 
for the D I profile contained in the data. 
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Fig. 10.— Fe II A2344 and A2382 line profiles from HST / STIS E230H spectra of BD +28° 4211 are compared 
with the best fit model and shown on a heliocentric velocity scale. The smooth solid line is the best 
fit obtained with three components. The histogram line is the normalized spectrum. 70% of the total 
column density is in the dominant component (right), 27% in the middle component, and 3% in the weakest 
component (left). The two stronger components are separeted by less than 8 km s _1 and are not resolved 
in Fe II features present in the FUSE spectra of BD +28° 4211. 
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Velocity (km s l ) 

Fig. 11. — Profile fits to three optically thin O I lines (AA919.92, 925.45, and 930.26). The top panels shows 
the fits for the SiC2 channel LWRS aperture spectra. The bottom panel shows the same transitions for the 
same channel but MDRS aperture. The definitions of the different line styles are the same as Figure 9. Both 
data sets are very similar, with a LSF FW HM ~ 20 km s -1 for LWRS and ~ 18 km s 1 for MDRS. 
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Fig. 12. — Profile fits to the N ] A952 multiplet. The top panels shows the fits to the spectrum for the 
SiC2 channel and the LWRS aperture. The bottom panel is the same channel for the MDRS aperture. The 
definitions of the different line styles are the same as Figure 9. The spectral resolution in the two data sets 
are very similar. Note the ER Lid — 0 P(3) line on the blue side of the triplet. This N I triplet is the only 
one available in the FUSE bandpass that is not optically thick and free of significant interference by lines of 
other species. 
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Table 1 

FUSE Observations of BD +28° 4211 


Program ID 

Date 

Aperture 

FP-split 
pos No. 

Offset 

(pixels) 

AA 

(A) 

Exposures 

Exp. Time 
( s ) 

M1080901 

Jun 13 2000 

LWRS 




1-4 

2192 

M1040101 

Jul 16 2000 

LWRS 

1 

-27 

-0.17 

1-14 

6251 




2 

0 

0.00 

14-25 

5073 




3 

+7 

+0.04 

26-30 

2454 




4 

+36 

+0.23 

31-36 

2799 

M1040105 

Sep 19 2000 

LWRS 

3 

+7 

+0.04 

1-8 

3691 




4 

+36 

+0.23 

9-17 

4146 

M1040102 

Jul 17 2000 

MDRS 

1 

-12 

-0.08 

1-13 

6279 




2 

0 

0.00 

14-30 

8211 




3 

+5 

+0.03 

31-42 

5795 




4 

+ 12 

+0.08 

43-52 

4540 


Table 2 

Properties of BD +28° 4211 


Quantity 

Value 

Reference 

Spectral Type 

sdO 

1 

l 

81.87° 

2 

b 

-19.29° 

2 

<F (pc) 

104 ± 18 

2 

V 

10.53 

3 

i U-B 

-1.26 

3 

B-V 

-0.34 

3 

Vhel 

+20.2 + 0.8 km s _1 

4 

T eff (K) 

82,000 

1 

log g (cm s~ 2 ) 

6.2 

1 

log(He/H) 

-1.0 

1 


a Trigonometric parallax 


References. — (1) Napiwotzki 1993; (2) Perry- 
man et al. 1997; (3) Kidder et al. 1991; (4) This 
paper (§4). 
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Table 3 

H 2 Column Densities' 1 


J 

No. Lines 

log(7V(H 2 ) ) 

0 

4 

14 45 

1 

8 

14 76 +0 16 
14 i / d - 0.06 

2 

13 

14 3Q+ 0 * 07 

3 

14 

14 17+ 0 * 05 
1 ( - 0.03 

4 

3 

< 13.68 


“Errors are 2<r. Value for J — 4 
is a 3 ct upper limit. 


Table 4 

D I and 0 I Equivalent Widths 3 


Line 

A 

(A) 

log /A 

W x (LWRS) 
(mA) 

W x (MDRS) 
(mA) 

(mA) 


D I 


SiCl 

SiC2 

SiCl 

SiC2 


Lyi 

920.713 

0.17 


10.4 ±4.7 


10.5 ±4.4 

10.5 ±3.2 

Ly/t 

922.899 

0.31 

11.1 ± 6.0 

13.7 ±3.8 

10.9 ±3.0 


12.3 ±2.2 

Ly»? 

925.974 

0.47 

18.5 ± 7.5 

18.3 ±7.0 

19.8 ±5.0 

16.0 ±4.0 

17.7 ± 2.7 

Lye 

937.548 

0.86 

34.4 ± 5.0 

30.3 ±5.7 

34.1 ±4.0 

29.1 ±5.0 

32.3 ± 2.4 

Ly<5 

949.485 

1.12 

42.0 ±4.6 

43.8 ±4.4 

48.5 ± 5.0 

44.0 ± 5.5 

44.4 ± 2.4 

Ly7 

972.272 

1.45 

50.6 ±4.0 

55.0 ± 4.4 

52.3 ±3.0 

54.6 ±5.0 

52.8 ± 1.9 


OI 


SiCl 

SiC2 

SiCl 

SiC2 



919.917 b 

-0.79 

15.6 ± 1.2 

20.8 ± 1.2 

16.7 ± 1.3 

16.6 ± 1.3 

17.5 ±0.6 


925.446 b 

-0.49 

27.4 ± 1.3 

23.7 ± 1.3 

28.1 ± 1.4 

23.1 ± 1.3 

25.7 ±0.7 


930.257 b 

-0.30 

32.5 ± 1.4 

34.0 ±1.4 

34.8 ± 1.5 

33.1 ± 1.4 

33.6 ±0.7 


919.658 c 

-0.06 

34.5 ± 1.3 

41.5 ±1.4 

43.1 ± 1.4 

30.9 ± 1.3 

37.2 ±0.7 


921.857 c 

0.04 

41.7 dh 1.4 

37.6 ±1.4 

41.3 ± 1.5 

33.4 ± 1.4 

38.3 ±0.7 


971.738 c 

1.13 

64.5 ± 1.5 

65.3 ± 1.5 

68.6 ± 1.6 

63.6 ± 1.5 

65.4 ± 0.8 


“Errors in this table are ]<r. 

b Singlet. Used /-values from Morton (1991). 

c Triplet. Used sum of /values from Morton (2000, private communication). 
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Table 5 

N(0 I) from Curve of Growth 1 


Channel 

Aperture 

log(JV(OI)) 
( cm 2 ) 

b 

( km s -1 ) 

SiCl 

LWRS 

16.22 ±0.10 

4.46 

SiCl 

MDRS 

16.25 ± 0.06 

4.83 

SiC2 

LWRS 

16.29 ±0.16 

4.41 

SiC2 

MDRS 

16.18 ±0.18 

4.24 


a Errors in this table are 2<r. 


Table 6 

D I , N I, and O I Profile Fitting Results' 1 


Case 

No. 

Windows 

DOF 

x 2 • 

Amin 

b 

( km s -1 ) 

log 7V(D I) 
( cm- 2 ) 

log N { N I) 
( cm" 2 ) 

log N { O I) 
( cm -2 ) 

Case l b 

39 

1444 

1674 

4.0 

14.99 ± 0.05 

15.51 ±0.05 

16.22 ±0.09 

Case 2 C 

39 

1405 

1441 

4.2 

14.99 ±0.07 

15.57 ±0.09 

16.22 ±0.16 

Case 3 d 

41 

1527 

1929 

3.8 + 3.3 

14.96 ±0.05 

15.49 ± 0.05 

16.17 ±0.06 


a Errors in this table are 2a. 

b One Gaussian component, LSF fixed for each channel (SiCl or SiC2). 
c One Gaussian component, LSF free to vary from window to window. 

d Three Gaussian components as defined by Fe II velocity structure, LSF fixed for each 
channel (SiCl or SiC2). See text. 


Table 7 

D I , N I, and O I Column Densities a 


Ion 

Profile Fitting 

Mean COG 

Adopted 

Solution 

log(7V(DI)) (cm" 2 ) 

14.99 ±0.03 

14.99 ±0.10 

14.99 ±0.05 

log(iV(N I) ) ( cm 2 ) 

15.55 ±0.13 


15.55 ±0.13 

log(7V(0 I)) ( cm -2 ) 

16.21 ±0.10 

16.23 ±0.08 

16.22 ±0.10 


Errors in this table are 2a. 


Table 8 

Interstellar Abundance Toward BD±28° 4211 a 


Quantity 

Value 

log N{D I) 

14.99 ± 0.05 cm- 2 

log 7V(N I) 

15.55 ±0.13 cm -2 

log N ( 0 I) 

16.22 ±0.10 cm- 2 

log AT(H I) 

19.8461^034 cm 2 

D/H 

1.39 ± 0.21 x 10" 5 

N/H 

5.08 ± 1.66 x 10“ 5 

O/H 

2.37 ±0.55 x 10" 4 

D/N 

0.274 ± 0.092 

D/O 

0.0590 ± 0.0133 


Errors in this table are 2<r. 




